Full counting statistics of super-Poissonian shot noise in multi-level quantum dots 
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We examine the full counting statistics of quantum dots, which display super-Poissonian shot 
noise. By an extension to a generic situation with many excited states we identify the underlying 
transport process. The statistics is a sum of independent Poissonian processes of bunches of different 
sizes, which leads to the enhanced noise. The obtained results could be useful to determine transport 
characteristics in molecules and large quantum dots, since the noise (and higher cumulants) allow 
to identify the internal level structure, which is not visible in the average current. 



Current fluctuations in mesoscopic structure far from 
equilibrium provide a great deal of insight into the rel- 
evant transport mechanism (see Refsi and£ for recent 
reviews). A complete understanding, however, requires 
to go beyond the second-order correlation function ('shot 
noise') and to study the full counting statistics 3 -, which 
yields all zero-frequency current-correlation functions at 
once. In contrast to transport of classical particles, which 
leads to Poisson statistics of the transfered charge, quan- 
tum transport of electrons is described by a binomial 
statistics 3 , which explains the suppression of noise in 
non-interacting conductors^. One way to observe super- 
Poissonian noise is to use a superconducting injector, 
which yields an enhanced noised as result of the doubled 
charge transfer due to Andreev reflection^. In systems 
with two superconducting leads even larger charge trans- 
fers can occust On the other hand, repulsive interactions 
as for example encountered in simple quantum dots can 
lead to a suppressed shot nois o?'?' 2 ? , similar to noninter- 
acting double tunnel junction systems. 

Some recent works found, surprisingly, a strong en- 
hancement of the shot noise in the transport through 
interacting quantum dots, in which the spin degeneracy 
was liftediS*ii. A similar enhancement was found in sys- 
tems of multiple coupled quantum dotsA 2 , in resonant 
tunneling diodes^ 3 - or due to cotunneling— . In addition, 
in multi-terminal quantum dots the enhanced noise re- 
sults in positive cross correlations^, which is probably 
the simplest mechanism to produce unusual positive cross 
correlations in a fermionic circuit^. We note here that 
these mechanisms do not require a negative differential 
resistance, which has been discussed in a number of works 
as possible source of enhanced shot noised. 

In this article we examine the full counting statistics of 
transport through a quantum dot, which contains several 
levels. We limit the discussion to the sequential tunnel- 
ing regime and for the moment to two energy levels e±, 
which in equilibrium are below the Fermi energy of the 
leads. Under the condition of strong Coulomb blockade, 
viz. the charging energy E c 3> e + — e— , the occupation of 
a single level may block the transport through the other 
level. This mechanism was termed dynamical channel 
blockade and can lead to a strongly enhanced nois o 11 ! 17 
(related mechanisms, due to coupling to phononsA& or 
bistable systems^, were recently discussed). Below we 



will show that the strongly enhanced Fano factor is a 
result of a combination of many different Poissonian pro- 
cesses, which transfer multiple charges. It follows from 
an analysis of the full counting statistics that these trans- 
port properties can only be understood from a study of 
all cumulants of the current. These transport proper- 
ties are markedly different from those in the situation, in 
which the levels are situated above the Fermi energy in 
equilibrium. Here, the fluctuations are similar to the non- 
interacting case and the noise is always sub-Poissonian. 

To study the full counting statistics, we recall here 
that we are dealing with the probability P{N) that N 
charges are transfered in a period to, which is related to 
the cumulant generating function S(x) as 
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The knowledge of the CGF is equivalent to the knowledge 
of the counting statistics. From the CGF we can obtain 
cumulants C n = —(—id x ) n S(x)\ x= o, which we identify 
with the current I = eCi/t , the noise Si = 2e 2 G^/to 
and similar for higher cumulants. 

The simplest system, which exhibits the phenomenon, 
is a two-terminal quantum dot with more than one level 
in the situation shown in Fig. QJl). In fact, this is the 
situation, which was first studied inA£. The applied bias 
voltage and the level configuration is such that the level 
e_ is located below chemical potentials of both leads and 
the level e+ in between. The charging energy is so large 
that only one of the levels can be occupied at the same 
time. 

In this limit the system is described by a Master equa- 
tion of the form 
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Here P± describes the occupation of the upper/lower 
level and Pq is the probability that the island is empty. 
The rates are parameterized by the bare tunneling rates 
Fl/r t° tunnel out of the dot to the left or to tunnel 
onto the dot from the right. The (small) parameter x is 
determined by the thermal occupations in the left lead: 
tunneling from the lower state is suppressed by a factor 
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FIG. 1: a) Two-terminal quantum dot with two levels. The 
tunneling rate from the lower level to the left lead is ther- 
mally activated due to finite occupation of the left lead. All 
other tunneling processes have much larger rates. Coulomb 
charging allows only one of the levels to be occupied, b) The 
same with N+ upper levels and iV_ lower levels, c) Trans- 
port cycles leading to super-Poissonian shot noise. With a 
small tunneling rate xYl a cycle is initiated by an electron 
tunneling from the dot to the left lead, followed by n cycles 
(n = 0, 1, 2, . . .) of tunneling through the upper level with 
probability 1/2. Finally the process is stopped by a tunneling 
event from the right terminal to the lower level with probabil- 
ity 1/2. In total n + 1 charges are transfered in a cycle with a 
rate a;r_L(l/2)" +1 , which explains the super-Poissonian noise 
and all other transport characteristics. 

x = 1 — /(e_), since the states in the left lead are mostly 
occupied. 

Now we determine the counting statistics using the 
method of Bagrets and NazarcA We add the corre- 
sponding counting factors exp(«x) to the entries of the 
Master equation, which involve changes of the dot occu- 
pation and refer to tunneling into one of the terminals. 
This procedure leads to the rate matrix 

/ ~xY L Y R \ 
M( X ) = -r L Y R , (3) 
\xY L e^ Y L e^ -2Y R J 

where we have chosen to count the charges in the left 
terminal. The counting statistics is determined from the 
lowest eigenvalue Ao(x) of M(x) according to 

S{x) = -Ao(x)io • (4) 

The eigenvalue equation reads 

= (xT L +X)[T L T R e^-(T L + X)(2T R + X)] 

+xT L e^T R (T L + A) . (5) 

For the present purpose, it is sufficient to determine the 
lowest eigenvalue perturbatively in x <C 1 and to assume 
Tl "C T^. The counting statistics is obtained as 

S( x ) = -2xr L t ^-±. (6) 



This cumulant generating function reproduces, of course, 
correctly the first cumulant C\ = 2xT '^to and the noise 
C*2 = 3Ci, corresponding to a Fano factor F = C2/IC1 of 
3. 17 The first interesting observation is that the current 
is twice the estimate of the simple thermal tunneling rate 
xTl- Furthermore, the noise corresponds to an effective 
charge q g = C%jC\ — 3, which could be interpreted as a 
Poissonian process of 3 charges. However, we obtain the 
for the higher cumulants 

C 3 = 13d , C 4 = 75Ci, (7) 
C 5 = 541Ci , C 6 = 4683Ci . (8) 

Notably, third (and higher) cumulants obey the relation 
Cn > q"g 1 Ci and, consequently, the transport statis- 
tics can not be explained by a simple Poissonian pro- 
cess in which multiple charges of size g e ff are transfered. 
Note that this would follow from a counting statistics 

-(xT L I g G ft)(exp(zg e ffx) - !)• 

A simple physical picture emerges, if we expand the 
cumulant generating function 10 in terms of Poissonian 
processes. As a result we find 

00 1 

s( x ) = - x r L J2^i^-i) . (9) 
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The counting statistics is therefore a sum of independent 
Poisson processes, in each of which a charge of ne is trans- 
fered as signaled by the factor exp(inx) — 1- Each process 
is weighted with a probability (1/2)". This result for the 
statistics suggests the following interpretation (see also 
Fig. The lower level is occupied most of the time, 
since it is well below the chemical potentials of the two 
leads. Due to Coulomb blockade the other level cannot 
be occupied at the same time (we assume here that the 
charging energy is larger than the bias voltage). At a 
finite temperature, there will be a small rate ~ xTl for 
the electron to hop out to the left lead. Next, the dot will 
be filled again with an electron from the right lead. How- 
ever, as we are in the situation eV R > e+ both levels will 
be occupied in the next step with the same probabilities 
of 1/2. In case the lower state is occupied we are back 
to the initial state and one charge has been transfered in 
this cycle, which is exactly what is happening in a single 
level quantum dot in the thermally activated transport 
regime. However, if the electron has tunneled onto the 
upper level, it can quickly tunnel out to the left lead and 
the dot is empty again and another electron can tunnel 
from the right terminal. Each intermediate cycle occurs 
with a probability 1/2 (since two states are available) and 
transfers an additional charge. Therefore each full cycle 
has a probability (1/2)" and transfers ne charge. 

We now turn to the more general situation of a multi- 
level quantum dot. The level structure we assume here 
is depicted in Fig.^ 3 ). As before the Coulomb energy is 
assumed to be so large that only one of the levels can be 
occupied at the same time. The levels are bunched into 
two groups, the lower group with iV_ levels and the upper 
group with N+ levels. We will furthermore assume that 
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the tunneling rate into and from all levels inside a group 
are the same. Physically this means that the maximal 
level spacing inside the lower group is smaller than the 
temperature, so that the thermal blocking factor x — 
1 — /(e~) is the same for all levels. For the upper group 
the requirement on the level spacing is less restrictive. 
We only have to assume that the band width of the upper 
group is smaller than the bias voltage. 

In principle we have to account now for N + + N\ + 1 
probabilities in the Master equations, but due to the 
simplifications made above, we can replace blocks con- 
taining 'upper'('lower') probabilities by a single equation 
with the rate multiplied with a factor A^+r-) , respectively. 
This procedure leads to the x-dependent rate matrix 



-xT L N- T R N+ 
M{ X ) = ( -T L N + T R N + 

xT L N-e l x T L N+e lx -2T R N+ 



(10) 



by n — 1 fast cycles, in which an electron tunnels through 
the dot with probability p. Finally the fast cycle stops 
with a probability 1 — p. As can be seen from Eq. i|12|) 
processes of all orders contribute to the transport char- 
acteristics. A simplified description in terms of a Pois- 
sonian process, in which an effective charge tunnels, can 
never reproduce the observed cumulants. It is worth to 
emphasize, that only the complete determination of all 
cumulants (or equivalently the full counting statistics) 
can unambiguously identify this transport mechanism. 

At this stage, we would like to remark that the sit- 
uation, discussed here, can only occur if the levels for 
zero bias are situated below the Fermi level (assuming 
symmetric capacitive coupling). In case the levels are 
initially above the Fermi level, the current and the noise 
show the usual behaviors 20 . The corresponding counting 
statistics can be obtained using the method of Ref. and 
turns out to be 



As before, it is sufficient to evaluate the CGF to lowest 
order in x. Let us introduce p = N+/(N- + N+), which 
is the relative probability that an electron tunneling from 
the right electrode occupies the upper level. The count- 
ing statistics (in first order in x) is 



S(x) = ^o(x)io = -~xT L t 
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This CGF gives the first cumulant C\ — xTLto/(l — p) 
and the higher cumulants are 
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We note that the average is larger by the factor 1/(1 —p) 
than the Poissonian estimate. Furthermore, higher cu- 
mulants obey the relation C n > (C2/Ci) ra-1 Cij2i which 
shows the unusual characteristics of the transport. Simi- 
lar to the previous case, we can decompose the FCS into 
multiple Poisson processes as 

oo 

S{x) = -xT L {l -p)^^ 1 (e mx - 1) • (12) 

71=1 

The statistics is a sum of Poissonian processes of multiple 
charges, weighted with probabilities P n = p n ~ 1 {l — p). 
This suggests the picture depicted in Fig. [fl;). The pos- 
sible Poissonian processes result from a tunneling event 
with the small initial rate xT^. This process is followed 



S(X) = g Vl +PW* - l)+PR(e-^ - 1) - ■ (I 3 ) 

Here pl,r are the respective probabilities of a tunneling 
event to the left (right) and T refers to the total tunneling 
rate. In the regime of unidirectional transport pi > p R 
and we obtain 



Ci=t -PL and F = l-p L /2. 



(14) 



The Fano factor is always between 1/2 for the symmetric 
case (pl — 1) and 1 for the asymmetric case (pl <1). It 
follows also, that higher-order cumulants arc suppressed 
even further. This situation is remarkably similar to 
transport through a noninteracting resonant level in the 
limit of a bias voltage much smaller than the width of 
the resonance. 

The difference between super-Poissonian and sub- 
Poissonian statistics discussed here has several interest- 
ing consequences. First, we note that the second cumu- 
lant in the super-Poissonian case is enhanced by a factor 
1 + 2N+/N-. If for example N- = 1, the enhancement 
of the noise is a direct measure of the number of excited 
states within the transport window. Notably, the full in- 
formation requires to study also the third (and higher) 
cumulants, since higher cumulants are not related to the 
second cumulant in a simple way. These observations are 
important for experimental studies of transport through 
molecules and other quantum dots with an unknown in- 
ternal level structures. As an example we consider a 
molecule, in which it is not known, whether transport 
takes place through the lowest unoccupied molecular or- 
bital (LUMO) or the highest unoccupied molecular or- 
bital (HOMO), because the position of the Fermi en- 
ergy of the leads is not known. The current in both 
cases shows a step at a threshold voltage. However, the 
noise is able to distinguish both situations. For trans- 
port through a multitude of LUMO levels, the noise is 
always sub-Poissonian, see e.gi&SS. However, transport 
through the HOMO falls in the category discussed in 
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this work. If there are excited levels in the transport 
window, which cannot be occupied in equilibrium due to 
Coulomb repulsion, we predict a super-Poissonian noise. 
Thus, noise measurements may be used to distinguish 
transport through the HOMO or the LUMO, which show 
the same average current- voltage characteristics. 

In summary, we have pointed out the unusual noise 
characteristics of the transport through a quantum dot 
with an internal level structure. A super-Poissonian noise 
results from a combination of Poissonian processes with 
multiple charges weighted with probabilities, depending 
on the internal level structure. The obtained results allow 
to identify the level structure, which is a useful informa- 



tion in molecular electronic devices. We emphasize here, 
that only a study of the complete statistics allows a de- 
termination of the transport mechanism underlying the 
super-Poissonian shot noise. While we have explored a 
rather simplistic model of a quantum dot, further stud- 
ies will reveal in more detail what information can be 
obtained from current fluctuation measurements in quan- 
tum dots and molecules with a complicated internal level 
structure. 
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That this relation is obeyed for all cumulants can be 
seen by writing C 2 (x) = Ci(x)q(x), where C n (x) = 
-(-id x ) n S{ X ) and q( X ) = (1 + pe' x )/(l - pe lx ) yields the 
effective charge q(0) — (1 + p)/(l — p). Since all deriva- 
tives w.r.t. ix of q(x) are positive, the chain rule yields 
Cn+i(0) > g(0)C„(0). 



